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A generalized Landauer formula, derived with the methods due to Keldysh, and Baym and 
Kadanoff, is gaining widespread use in the modeling of transport in a large number of different 
mesoscopic systems. We review some of the recent developments, including transport in semicon- 
ductor superlattices, calculation of noise, and nanoelectromechanical systems. 



I. INTRODUCTION 



■ A. Some background, and a few historical remarks 

A mesoscopic transport measurement is often concerned with a situation where the " device" , whose properties are 
the subject of the investigation, is connected to structureless "contacts" via "ideal leads", i.e., in a situation which is 
accessible via the Landauer formula (or, more generally, the scattering approach to transport). An important feature 
is the fact that the size of the device is finite, and comparable to other important length-scales of the system, such as 
the phase-braking length or impurity mean free path. Thus the wave-properties of the charge carriers are important, 
leading to a number of interesting interference effects, such as weak localization, or universal conductance fluctuations. 
1 1 The conductance g of the device, for example, is then given by the celebrated Landauer formula[l, 2] 
■ 

S=-H*l> (!) 

where t is the quantum mechanical transmission amplitude through the device. 

The expression given above holds for the one-channel case, for low applied voltages, and in a situation when the 
"device" can be modelled by a noninteracting system. It is then natural to ask: Can this equation be generalized 
— ■ to the case, when interactions are important? Or: Can this equation be extended to time-dependent situations, or 
situations where superconductivity or magnetism are essential for the physics? How about strong driving fields? 

Many authors have addressed these questions, with a number of different approaches. In the spirit of this conference, 
this review is concerned with the subset of these theories which use the nonequilibrium Green function technique. The 
earliest applications to mesoscopic transport, to my knowledge, are due to French researchers [3-6]: these researchers 



were mainly interested in inelastic effects in tunneling through oxide barriers. It is a curious side-note to observe that 
these early, and pioneering, papers were essentially forgotten during the 80 's, but have obtained a substantial revival 
■ since mid-90's, and are presently cited more often than ever earlier. The explanation lies perhaps in the fact that 
the whole idea of mesoscopics is newer than these early papers, and it took the mesoscopic community a few years to 
1 ■ realize the applicability of these ideas. 

For the purposes of the present review, the next important development was the paper by Meir and Wingreen[7], 
which gave a very useful formal expression for the current in terms of the exact Green function of the device (or 
"central region"). This formula, and similar expressions obtained by other groups, were then applied to the Kondo 
problem out of equilibrium, a notoriously difficult problem, which remains a topic of active research even today. My 
previous review[8] in the first meeting of the present series focused in some of these issues, paying particular attention 
to the time-dependent generalization of the Meir-Wingreen expression [9]. I shall not repeat any of that material, 
but rather focus on other topics: the selection criterion has been that either they were not discussed during the 
first meeting, or that they are strictly post-1999 vintage. I have chosen to discuss three examples: (i) Transport 
in a semiconductor superlattice; (ii) Calculation of the noise in a spintronic system; and (iii) Tunneling transport 
in a nanoelectromechanical (NEMS) device. These three topics have a common feature: they all have practical 
applications, and, dare I say in the present meeting, even commercial potential. 

It should be noted that there are many other recent applications of the NGF to mesoscopic transport which, due 
to space and time limitations, the present review does not address. One such example is transport in nanowires, 
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fabricated either with scanning tunnelling microscope, or with break junctions, which is a large research field where 
many theoretical calculations emply the NGF techniques. An exhaustive review has recently become available [10]. 
Molecular electronics holds enormous potential, and here a combination of ab initio electronic structure calculations 
(within the density- functional scheme) and NGF appears to the most promising theoretical technique [11, 12]. Yet 
another example is the "circuit theory" developed by Nazarov[13, 14], which has successfully been applied to a number 
of hybrid structures, consisting of superconductors, ferromagncts, or semiconductors. The more abstract field-theoretic 
formulations, based on path-integrals and/or Grassmann variables fall also outside our present purposes, even though 
they play an important role in the study of disordered systems, or dephasing due to the environment. 



B. The basic equations, and their limitations 



For completeness, we sketch here a derivation of the basic expressions used in the theory. Several more complete 
accounts are available elsewhere[8, 9, 15]. A brief reminder of how the nonequilibrium formalism works in the 
context of mesoscopic transport measurements is also in place. One reasons as follows. In the remote past the 
contacts and the central region (i.e., the "device") are assumed to be decoupled, and each region is in thermal 
equilibrium. The equilibrium distribution functions for the three regions are characterized by their respective chemical 
potentials; these do not have to coincide nor are the differences between the chemical potentials necessarily small. The 
couplings between the different regions are then established and treated as perturbations via the standard techniques 
of perturbation theory. The nonequilibrium nature of the problem manifests itself in that symmetry of remote past 
and remote future has been broken, and thus one must do the calculations on the two-branch time contour. It is 
important to notice that the couplings do not have to be small, e.g., with respect level to spacings or fceT, and 
typically must be treated to all orders. 

Let us next consider some generic Hamiltonians: H = Hl + Hr + Ht + H CCD , or, explicitly: 



k,a£L/R 

+ X] [ V ka:ncl a d n + Il.C. 
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+H ccn [{d n },{di}} , 



(2) 



where the central part Hamiltonian must be chosen according to the system under consideration. The operators 
{d n }, {d^} refer to a complete set of single-particle states of the central region. Occasionally we specify explicitly 
the orbital and spin quantum numbers: n = m,a, and analogously for the states in the leads. The derivation of the 
basic formula for the current does not require an explicit form for H cen ; the actual evaluation of the formula of course 
requires this information. We write H ccn = J2 n e„d^c?„ -I- H lnt , where H lnt could be electron-phonon interaction, 
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or an Anderson impurity: 

H L = U Yj d ™,1 dm '1 d ln,l d ™>l ■ 

m 

The current operator for the (say) left lead is 



(3) 
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The physically relevant observables can be expressed in terms of expectation values of the current operator, or its 
higher powers. For example, one can show[9, 15] that the current leaving the left contact is 



(I L ) = J L {t) = ~J* dtxj ^ImTr{e- i <* 1 -*)r I '(e,i 1 ,i) 
x [&<(*,*!) + /£(e)G'"(t, ti)] } . 



(6) 
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Here the Green functions are defined by 



G< m (t,h) = iidl^dnit)) (7) 

G r nm (t,h) = -i0(t-t 1 )<[dn(*),4(*i)]>. (8) 



T m „ describes the coupling between the central region and the contacts, and /£(e) is the equilibrium distribution 
function of the left contact. In the dc-limit, (6) reduces to the result of Meir and Wingreen[7]: 

+ [f° L (e)T L (e) f R (e)T R (e)} [G r (e) G»(e)j } (9) 

= J J ^[h{e)-f R {e)]T{e), ( 10 ) 

where 

T(£) = Tr { r'(ffr^) [ Gr '-''- G- ^}- <"» 

The expressions (6) and (10) are the central formal results whose consequences we explore in this review. [53] They are 
formally exact, and give the tunneling current for an interacting system coupled to noninteracting contacts (or, more 
precisely, for contacts which can be described by an effective single-body Hamiltonian). Thus, in a time-dependent 
situation the displacement current must be considered separately. It should also be noted that these equations only 
define the starting point of any calculation: to get into physical results one must evaluate the correlation function 
and the retarded/advanced Green function, Eqs.(7) and (8), respectively. These functions obey the Kcldysh equation, 
and the (nonequilibrium) Dyson equation: 

G< - G r S<G a , (12) 
G r = G r + G r Z r G r . (13) 

The success of the theory depends on whether one can construct a self-energy functional that captures the essential 
physics, and that a good solution can be found for these coupled equations. Both of these steps may be hard indeed. 



II. TRANSPORT IN A SEMICONDUCTOR SUPERLATTICE 



In 1970 Esaki and Tsu[16] suggested that semiconductor superlattices, man-made structures which consist of alter- 
nating layers of different semiconductor materials, would have physical properties which could be used for a number 
of device applications. Very shortly, the spatial variations in the band-gaps will lead to a spatially varying conduc- 
tion band edge, which supports minibands, which in turn diplay very interesting transport properties, such as Bloch 
oscillations, or negative differential resistance. A well known-result is the Esaki-Tsu IV-characteristic, 

i(v) = 2/ max v ^qr^> ( 14 ) 

where I max and Vq depend on the physical paratmeres of the system, such as the superlattice period, scattering 
rate, and temperature. To derive expressions like this, three main approaches have been used in the literature: 
(i) Miniband transport[16], (ii) Wannier-Stark hopping[17], and (hi) sequential tunneling [18]. The three different 
approaches have different domains of validity, and are all likely to fail if the three basic energy scales, i.e. scattering 
induced broadening, miniband width, and potential drop per period all have comparable values. The basic features 
of these three approaches are summarized in Figure 1. 



In order to map out the boundaries of the various domains of validity, and to access the region where the approaches 
(i)-(iii) fail, a higher level theory is required. To achieve this, Andreas Wacker and myself, together with several 
colleagues, launched a program whose task was to develop a nonequilibrium Green function theory for superlattice 
transport. [54] Certain aspects of this program are now complcted[22, 23], and in what follows I will review some of 
the highlights. It should be noted that the literature on superlattice transport is vast and here I can give only a 
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TABLE I: The three standard approaches to miniband transport, and the physical picture underlying them. (Courtesy of A. 
Wacker.) 



Miniband conduction 



Wannier Stark 
hopping 



Sequential tunneling 



coupling Ti 



exact 
miniband 



field drop eFd 



acceleration 



exact: Wannier Stark states 



lowest order 



energy 
mismatch 



scattering T 



golden rule 



golden rule 



"exact" 
spectral 
function 



very superficial discussion; the reader is referred to two recent review articles where a much fuller account can be 
found[24, 25]. 

Let me start with a few disclaimers. The quantum theory has not yet been fully developed to the case when the 
electric field is inhomogeneous (domain formation), nor is it available for the time-dependent case (photo-assisted 
transport; progress is however being made see, e.g., Appendix C in the review by Wacker[24]). For these important 
situations one has to apply one of the simpler approaches discussed above. As far scattering is concerned, impurity 
scattering and phonon scattering have been discussed, but carrier-carrier interaction is still a future task. 

The task is now to solve the coupled Keldysh and Dyson equations, Eq. (12-13). We adopt the tight-binding 
representation of the single-particle Hamiltonian: 



i)Ti + S n>m (E k - neFd), 



(15) 



where T\ is the nearest neighbor coupling, Ek = H 2 k 2 /(2m) the kinetic energy perpendicular to the growth direction, 
F the applied field, and d the superlattice period. In this basis the Keldysh and Dyson equations (12-13) read 



Gfnn(E) = J2 G ^ [E + eFd 



E + eFd mi 



,mi — n 
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(16) 
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(17) 



Next one needs to specify the self-energies. We have considered [2 2, 23] impurity scattering, optical phonon scatter- 
ing, and mimicked acoustic phonon scattering by a very low-energy optical phonon, all in the self-consistent Born 
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FIG. 1: The range of validity of various approaches to superlattice transport, in the parameter space spanned by the nearest- 
neighbor coupling Ti, and the potential energy drop edF per period, in units of the scattering width V. 



A=20.3 meV T=77 K 




10 4 10 3 10 2 10 1 
eFd [eV] 

FIG. 2: Drift velocity vs. applied field. 

approximation. By numerically solving these coupled equations, computing the current, and comparing to the cor- 
responding IV-curves found by the simpler approaches (i)-(iii), we can construct a "phase-diagram" (see Figure 2), 
which indicates where the simpler approaches hold, and where a quantum approach is necessary. 

We have also compared the quantum mechanical drift velocity vs. field relation to the results obtained with a 
semiclassical Monte Carlo simulation. This is quite interesting because the two methods are totally different, and 
both require computationally rather intensive calculations. Typical results are shown in Figure 3. For the parameters 
considered here, the Monte Carlo simulation gives very good results, except that it misses the weak phonon replica 
seen in the quantum calculation. 

The approach sketched here is ideally suited to transport phenomena where quantum phenomena, such as resonant 
tunneling, or phonon-assisted tunneling, play an important role. Another application concerns quantum cascade 
lasers[26], where the current injection occurs through a "funnel": the superlattice is designed so that the miniband 
width varies with distance. Another recent calculation concerns the evaluation of gain in such structures[27]. 

III. NOISE IN SPINTRONICS 

The emerging field of spintronics [28-30], where, in addition to the charge, also the electron spin is used to design 
new devices, has led to fascinating and novel ideas such as spin filters [3 1-33], spin field-effect transistors [34], and 
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proposals for solid state quantum computing[35]. For example, quantum dot systems can in principle be used to 
control the electron spin and are thus suitable for creating quantum bits relevant for quantum gate operations [36]. 

A detailed theoretical study of nonequilibrium transport properties of spintronic devices is necessary in order to 
understand the basic physical phenomena and to predict new functionalities. Calculation of the current, for example, 
can give the conductance/resistance of a system and its dependence on magnetic field, Coulomb interaction, spin-flip 
and so on. On the other hand, current fluctuations, due to the granularity of the charge (shot noisc[37]), are also 
relevant because their measurements can provide additional information not contained in the average current [38]. 

Here we illustrate how the nonequilibrium Green function technique can be used to calculate current and its 
fluctuations (noise) [55] in a quantum dot coupled to two ferromagnetic leads as a function of the applied voltage for 
parallel (P) and antiparallel (AP) lead-polarization alignments. We include Coulomb interaction in the Hartree-Fock 
approximation as well as spin- flip in the dot. We show that spin- flip makes the alignment of the lead polarizations 
less important; both P and AP results coincide for large enough spin-flip rates. This fact gives rise to a reduction of 
both Fano factor [56] and tunnelling magnetoresistance (TMR) as we show below. 

We model the central region with the Hamiltonian 

H D =^2 e o d U<? + Un ]ni + R(d\d l + d{d T ), (18) 

(7 

where d a (<4-) destroys (creates) an electron in the dot with spin o and the energy eo is spin independent [43, 44]. 
In addition, wc assume that the dot is a small enough in order to have only one active level eo- In the presence of 
a voltage the level shifts by eo = e^ — where e^ is the dot level for zero bias. In a more realistic calculation 
one should determine the bias dependence self-consistently. The spin-flip scattering amplitude R is viewed here as a 
phenomenological parameter. The spin-flip process lifts the degeneracy, splitting the quantum dot level to two states, 
let us call them ei^, with corresponding operators. The current is readily evaluated with the formulas given in Section 
1.2 with the result 

J L = yRe J dt 2 Tr{G r {t,t 2 )Z L< (t 2l t) + G<(t,t 2 )Z La (t 2 ,t)}, (19) 

where G r and G< arc the nonequilibrium dot Green functions, with elements Gfj(t, t 2 ) = i(dj(t 2 )di(t)) and G\j{t, t 2 ) = 
—i6(t — t 2 )({di(t), dj(t 2 )}). The lesser (retarded, advanced) tunnelling self-energy is given by 

Z L <™ (*2,i) = J^>Ll (20) 

z k 

x (9S' a \t 2 ,t)+ 9 <l' a \h,t) 9^' a) (t 2 ,t)- 9 <^\t 2 ,t)\ 

^? o) (*>.*)-^r o) (fa,*) d^M + a^MJ ' 

where g^a^ i s tne l esser (retarded, advanced) uncoupled Green function for lead L. Equation (20) leads to a 
generalization of the coupling T found in Section 1.2 above; the coupling matrix now becomes 

At + r| rf-rf\ 

~ 2 \ Tf -Tf Tf +Tf) ■ [21 > 

Accounting for Coulomb interaction in the Hartree-Fock approximation, we can write down a matrix Dyson equation 
for the retarded Green function, G r = G 0r + G 0r S r G r , and a Keldysh equation for the lesser Green function 
G< = G r S < G a , where G 0r is the uncoupled dot Green function. In these equations the self energies are the sum of 
the left and right self energies, i.e., £( r ' < ) = £- L ( r ><) + v^(»%<). A self-consistent calculation is required to calculate 
(nj) and (dldi), which are given by the lesser Green function, (djdj) = J |^ImG<(cj). 

The current operator can be written as its average value plus some fluctuation, i.e., I v = J v + SI V (here i] = L/R 
labels the contacts). In our system there are two sources of noise, namely, thermal noise and shot noise. The first 
one is due to thermal fluctuations in the occupations of the leads. It vanishes for zero temperature, but can be 
finite for T ^ and eV = 0. On the other hand, shot noise is due to the granularity of the electron charge; it is a 
nonequilibrium property of the system in the sense that it is nonzero only when there is a finite current (eV ^ 0). 
To calculate the noise (thermal+shot noise) we use the definition S m i (t — t') = ({SI n (t) , 8I v i (t')}), which can also be 
written as S m i{t — t') = ({I v (t), I v > (t 1 )}) — 2J^. After a lengthy but straightforward calculation [45], we find for the 
noise power spectrum (dc limit; a scalar version of this equation has been found earlier [42]) 
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(22) 



The dc noise (zero frequency) is position independent, and it is possible to show that Sxl(O) = Srr(0) — —Slr(0) = 
—Srl(0)[37]- In our numerics we make a few simplifying assumptions. We assume that the couplings T v are energy 
independent, but allow a polarization dependence. For the physical parameters we use accepted values from the current 
literature [43]. Our Hartrec-Fock approximation for the electron-electron interaction does not include correlations of 
the Kondo type, however we do not expect these to change our results in the present range of parameters. 



0.01 



— r— | i | i \^=f=n r 

(a) 




FIG. 3: Current and noise as a function of the bias for parallel (P) and antiparallel (AP) alignments and with R = and 0.1(7. 
The curves for R — 0.1U are only for the AP alignment; observe that these are almost on top of the P curves, except within 
the sloping region around U and 3U. Both current and noise are reduced when the right lead changes its polarization from P 
to AP, following the typical behavior of TMR. The inset shows a suppression of the AP-Fano factor due to spin-flip. 

Figure 4 shows current (a) and noise (b) as a function of the bias with R — (solid line) and R — 0.1U (dotted line) 
for both P and AP configurations. Because P and AP curves for R = 0.1U coincide, we plotted only the AP case. The 
first enhancement of the current and noise at eV — U happens when eo crosses the left chemical potential, allowing 
electrons to tunnel from the emitter (left lead) to the dot and then to the collector (right lead) . The current and noise 



8 



remain constant until the second level eo + U reaches /il at eV = 3U , when another enhancement is observed. Each 
enhancement corresponds to a peak in the differential conductance a^is- When the system changes from parallel (P) 
to antiparallel (AP) configurations the current is reduced. This is a typical behavior of tunnelling magnetoresistance 
(TMR). The noise is also affected by this resistance variation, showing a similar reduction. 

Looking at the effects of spin- flip on current and noise we see that the AP curves with R — 0.1U (dotted lines) tend 
to lie on the P curves with R — 0, thus showing that lead alignments are less important when spin-flip plays a part. 
This AP current enhancement due to spin- flip gives rise to a reduction of the TMR; since TMR = (Ip — Iap)/Iap, 
when Iap — ► Ip we have TMR — > 0. For a somewhat simpler model W. Rudzihski et aZ.[43] found a similar behavior. 

IV. NANOELECTROMECHANICAL SYSTEMS 

Microelectromechanical systems (MEMS) are today an important part of our technology. Their functionality is 
based on combining mechanical and electronic degrees of freedom, the great advantage being that the whole device can 
be fabricated with standard Si-processing technology. Typical applications include hearing aids, sensors, or actuators. 
As the fabrication technology gets refined, we soon expect to find systems where the mechanical parts are in the 
nanometer range, see the review by Craighead [47]; hence the acronym NEMS. An example is the nanomechanical 
electron shuttle constructed by Erbe et aZ. [48] (see Figure 5), based on the theoretical ideas of Gorelik et al.. [49] 
Transport through the device was modelled with rate-equations [48], however these are not expected to work when one 




FIG. 4: The nanomechanical electron shuttle. An ac-voltage coupled to the gates Gi and Gi causes the central electrode, the 
"clapper", G to oscillate, and charges tunnel from the source S to the island, and from the island to the drain D. 

enters the coherent transport regime, at lower temperatures and smaller devices, and a quantum theory of transport 
must be invoked. Such an approach has recently been formulated by Fedorets et al. [50] , and I'll give a brief introduction 
to this topic. 

The archetypal NEMS device consists of a moving part, and electrodes. We shall model the moving part by 
a quantum well, whose coupling to the electrodes are position dependent: the tunnelling amplitude is written as 
T L m(x(t)) = Tp/p exp[=px(t)/A], where A is some characteristic tunnelling length, and x(t) is the center-of-mass coor- 
dinate of the moving quantum dot. The electronic degrees of freedom are governed by exactly the same Hamiltonians 
as discussed in Section 2; now the dot-level depends on the dot's location via a = eo — Ex(t), where E is the electric 
field, which depends on the applied bias. [57] The center-of-mass of the dot obeys Newton's equation of motion, 

x + ujIx = F(x)/M , (23) 

where M is the mass of the grain, ujq — k/M is the characteristic frequency, and F is the force, which includes both 
the electric force acting on the charge(s) on the dot, and an exchange force, which arises due to the x-dependence of 
the tunnelling matrix elements. One can evaluate the force via F = —(dH/dx), and the result is 

F(t) = -iEG<(t,t) + 2/\J2(-irT a (t)lm[G< k (t,t)}, (24) 

a.k 

using the notation of Section 2. The problem is thus (again) reduced to the determination of the two lesser Green 
function: G < for the dot, and the non-diagonal G< k = i(c[ a d(t)). Fedorets et aZ. [50] calculate these for a nonin- 
teracting system, and proceed to present an analysis of the mechanical stability of the system: does the dot execute 
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regular oscillations (whose frequencies are determined by an appropriate linearization of (23) and (24), or does it 
perhaps become unstable, as the bias is increased? The details of the analysis are not our concern here; the upshot is 
that above a certain threshold value an instability results. The analysis has also bearing on a recent experiment [51], 
where vibronic anomalies were observed in a single-Ceo-transistor when current was passed through it. 

The analysis of Fedorets et al. [50] is very interesting and suggests for several further refinements. For example, 
what are the results for an interacting system (Coulomb blockade)? How about the environmental degrees of freedom? 
This issue was addressed recently Armour and MacKinnon [5 2], in a slightly different context. Finally, is it possible to 
combine the spintronic effects with charge shuttles? Can one envisage a spin shuttle? Will there be a new technology 
called NEMSS (nanoelectromechanical spin systems)? 

V. CONCLUSIONS 

I have reviewed some of the post-1999 developments in applying NGF to modelling of transport in mesoscopic 
systems. The common theme in my review has been to focus on "real devices", which may have "real applications". 
I find it very pleasing that the NGF technique, often regarded as an academic exercise most suited for theoretical 
games, is now becoming a strong tool in the analysis of practical devices. This trend is also confirmed by several other 
talks at this conference. At the same time there is still much room for theoretical refinements, and I'm convinced 
that in the coming years we will witness significant progress in this field, both abstract and practical. 
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